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Left-handed materials are studied numerically using an 
improved version of the transfer-matrix method. The trans- 
mission, reflection, the phase of the reflection and the absorp- 
tion are calculated and compared with experiments for both 
single split-ring resonators (SRR) with negative permeability 
and left-handed materials (LHMs), which have both the per- 
mittivity e and permeability fj, negative. Our results suggest 
ways of positively identifying materials that have both e and 
fi negative, from materials that have either fj, or e negative. 

PACS numbers: 73.20.Mf, 41.20.Jb, 42.70.Qs 



Photonic band gap (PBG) structures were originally 
introduced to control the electromagnetic, wave propa- 
gation of .HiatterEm. Not only dielectricElfl but metallic 
structurea3~Q were proposed for applications in the mi- 
crowave and infrared regions. Very recently, a new area 
of research, called left-handed-materials (LHM), has been 
introduced by Pendry et alol and Smith et al.@. LHM 
are by definition composites, whose properties are not de- 
termined by the fundamental physical properties of their 
constituents but by the shape and distribution of specific 
patterns included in them. Thus, for certain patterns and 
distribution, the measured effective permittivity e e ff and 
the effective permeability [i e g can be made to be less 
than zero. In such materials, the phase and group ve- 
locity of an electromagnetic wave propagate in opposite 
directions giving rise to a number of novel properties!! 
This behavior has been called "left-handedness" , a term 
first used by Veselagot3 over thirty years ago, to describe 
the fact that the electric field, magnetic intensity and 
propagation vector are related by a left-handed rule. 

By combining a 2D array of split-ring resonators 
(SRRs) with a 2D array of wires, Smith et al.cl demon- 
strated for the first tinae the existence of left-handed ma- 
terials. Pendry et al.13 has suggested that an array of 
SRRs give an effective /j, B s, which can be negative, close 
to its resonance frequency. It is also well knowrafl that 
an array of metallic wires behaves like a high-pass fil- 
ter, which means that the effective dielectric constant is 
negative at low frequencies. Very recently, Shelby et al.Eil 
demonstrated experimentally that the index of refraction 
n is negative for a LHM. Also, PendryEj has suggested 
that a LHM with negative n can make a perfect lens. 
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In this Letter, we present systematic numerical results 
for the transmission, reflection and absorption proper- 
ties of left-handed materials. An improved version of 
the transfer-matrix methodEj fEMM) is used. Qualita- 
tive agreement with experimentsEnlJ is obtained, and new 
predictions about the absorption and the phase of the re- 
flection, to be checked experimentally, are presented. 

The transfer-matrix methocO is used to calculate the 
EM transmission and reflection of LHMs. In the TMM, 
the total volume of the system is divided into small cells 
and fields in each cell are coupled to those in the neigh- 
boring cell. Then, the transfer matrix is defined by relat- 
ing the incident fields on one side of the LHM structure 
with the outgoing fields on the other side. The main ad- 
vantage of the TMM is the calculation of the transmis- 
sion and reflection coefficients for EM waves of various 
frequencies incident on a finite thickness slab of a LHM. 
In this case, the material is assumed to be periodic in the 
directions .parallel to the interfaces. The TMM has, been 
previouslya applied in studies of metallic structurestTElo, 
as well as in dielectric structures. In all these examples, 
the agreement between theoretical predictions and exper- 
imental measurements was very goodH .— . 

Following the algorithm described ir£3, we develop 
the new version of the transfer-matrix code. The main 
change is the faster normalization of the transmitted 
waves in the calculation of the transmission coefficient 
through the structure. For the structures described in 
Figure 1, the typical CPU time necessary to calculate 
transmission for a given frequency is, approximately 55 
sec/unit cell on an alpha workstationt^l Most of the CPU 
time is spent in the normalization of the actual numerical 
data. As our structure is highly non-homogeneous (per- 
mittivity of systems varies in three or more orders from 
site to site) , we have to normalize the numerical data af- 
ter 3-4 steps. For studies of more homogeneous samples, 
the CPU time could be reduced by a factor of 2 to 5. 

The TMM has been used to simulate the reflection 
and transmission from an array of square split-ring res- 
onators (SRRs). Figure la shows a diagram of a single 
square SRR|-af the type used for our simulations and for 
experimental^. In our case we have c = d = g = 0.33 
mm and the size of SRR w — 3 mm. Figure lb shows a 
three dimensional realization of the actual LHM that we 
have simulated. 

Figure 2 presents the results of the transmission versus 
frequency for the split ring resonators alone, and of the 
LHM, which consists of the SRRs with metallic wires 
placed uniformly between the SRRs. The square array 
of metal wires alone behaves as a high pass filter with a 
cutoff frequency v c = 19 GHz. The cutoff frequency v c 
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of the metallic wires is given by v c = cn g \ 1 t/2dy/eo, where 
d is the distance between the wires, ci; g ht is the velocity 
of light in the air and eo is the dielectric constant of 
the background. In our case d = 5 mm and eo = 1. The 
cutoff frequency is independent of either value of Re e m or 
Im e m of the metal, provided that either |Re e m | > 1000 
or Im e m > 1000. The dot-dashed curve is that of the 
SRR array with d = 5 mm. By adding wires uniformly 
between the SRRs, a pass band occurs where both /Li e ff 
and e e ff are negative (solid line). The transmitted power 
of the LHM is very high (close to one), because no Im 
e m is taken for the metal. We also studied more realistic 
cases with non-zero Im e m . We found that the resonance 
frequency z/ increases as |e m | increases and saturates to 
a value of vq w 8.5 as |e m | — > oo. This is clearly shown 
in Fig. 3 where we plot i>o versus the magnitude of e m = 
e r + iei. The resonance frequency is also sensitive to the 
permittivity of the board that the SRR is lying and of the 
embedding medium. We have found that vq drops as the 
dielectric constant of the embedding media e a increases. 
Vq drops from its value of 8 GHz for e a — 1 to a value of 6 
GHz when e a = 2.0. The value Vq increases by decreasing 
the value of the permittivity e;, of the dielectric board. i>q 
increases to a value of 9 GHz from 8 GHz when ef, drops 
to a value of 1.4 from its original value of tb = 3.4. 

The results presented in Fig. 2 are done by assuming 
that the permittivity of the metal e m had a large negative 
value. As one can see from Fig. 2, the value of the 
transmission for the LHM is very close to one. This value 
drops considerable if one uses the fact that e m for the 
metal has a complex permittivity. In Fig. 4a we present 
results similar to that in Fig. 2, with the only difference 
being the value of the Im e m . Notice that the position 
of the resonance peak is at the same value at around 7.8 
GHz. Due to the non-zero Im e m there is considerable 
drop of the peak of the transmission for the LHM. We 
have not done a systematic effort to fit the experimental 
results of Smith et al.a. 

In Fig. 4b, we present new results of the frequency 
dependence of the absorption for either the SRR or the 
LHM. For the arrays of SRRs alone, we find that absorp- 
tion has two peaks, shown as dashed lines in Fig. 4b, 
while the LHM has only one absorption peak (shown as 
a solid line in Fig. 4b) at the center of its resonance fre- 
quency. This might be a reasonable way of identifying 
materials that are left-handed and have both their e e s 
and /Lteff negative. On the inset of Fig. 4b we present 
the frequency dependence of the absorption as they were 
found for a simple model, usingEn3 the following forms 
of the frequency dependence of the effective permittivity 
and permeability: 
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where v p is the plasma frequency or cut-off frequency v c 
of the wires, v c = v p 19 GHz in our case. A negative 
/i e ff accounts .for the deep in the transmission in SRRs 
and its formBfl is given by 



where F is the filling factor which is in our case close to 
0.3 and vq w 7.9 GHz is the resonance frequency of the 
SRR. We used that 7 = V = 0.2 GHz. There is very 
good agreement between the solution of present model 
and numerical simulations of real structures. 

In Fig. 5 we present the results of the phase of the re- 
flection for both the LHM and the array of SRRs. Notice 
that there is a substantial difference in their character- 
istics, which can help us identifying meta-materials that 
either have negative /Lt e ff or both negative /j, e s and e e ff. 

Finally, in Fig. 6 we present results of our numerical 
experiments with various vacancies in SRR and LHME3. 
First, we remove one SRR from the structure. In PBG 
materials, where a gap is created by the interference effect 
of propagating waves, such a vacancy causes the appear- 
ance of sharp resonance frequency inside a gap. As it is 
seen in Fig 6a, the absence of one split ring resonator does 
not influences neither the position of the gap nor its form. 
Only the amplitude of the transmission decreases. This 
proves that the origin of the gap is not due to the interfer- 
ence effect but, as we assumed, in the negative effective 
permeability of the structure. Fig. 6 also shows how the 
transmission peak changes when one SRR or one wire is 
removed from the periodic structure of a LHM. In the 
absence of one SRR, the transmission peak decreases by 
more than one order of magnitude and becomes narrower 
(see Fig. 6a), while the wire vacancy creates an "inverse 
resonance" inside the transmission peak (see Fig. 6b). 

We used an improved transfer-matrix method for nu- 
merical studies of complex meta-materials. Our numeri- 
cal data of transmission confirms the presence of a reso- 
nant gap, in agreement with theoretical predictions and 
experiments. Our method enables us to analyze also the 
reflection and absorption of the light and also the change 
of the phase of the reflected field. Obtained numerical 
data are consistent with a simple homogeneous model 
defined by effective permittivity and permeability. This 
supports our belief that for the frequencies in the reso- 
nance gap our structure possesses negative effective /i e s 
and negative effective e e ff. Although our method does 
not allow us to vary the size parameters of SRR contin- 
uously, we are able to predict, at least qualitatively, how 
the position of resonance gap depends on various param- 
eters of the system. We believe that generalization of our 
numerical method to non-homogeneous discretization of 
space will enable us to analyze more realistic structures. 
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FIG. 1. (a) SRR consists from two splitted metallic rings 
located on a dielectric board. We approximate the rings by 
squares of size w. Other parameters, which influence the res- 
onance frequency vq are the ring width c, the radial gap d, 
and the azimuthal gap g. (b) Configuration of the unit cell of 
a LHM. Light propagates along the z axis and the magnetic 
intensity is parallel to the SRR axis, which is the z-axis. We 
assume periodic boundary conditions in x and y directions. 
In the numerical simulations the unit cell is 5 x 3.63 x 5 mm 
and is divided to N x X N y X N z = 15 X 11 X 15 small cells. 
This defines the length unit =0.33 mm. Number of unit cells 
in z direction varies from 6 to 20, and is 9 for most of the 
cases presented here. 



3 



6 7 8 9 

Frequency [GHz] 

FIG. 2. Transmission of electromagnetic waves through 
arrays of SRRs, wires, and through a LHM with a real per- 
mittivity for the metal e m = —3000. The dielectric constant 
of the board is ei = 3.4. The thickness of wire is 1 x 1 mm. 
Magnetic intensity is parallel to the SRR axis. For the other 
polarization, no resonant gap was observed. 
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FIG. 3. Dependence of the resonance frequency vq on the 
absolute value of the permittivity |e m | of metallic compounds. 
Both real (open circles) and complex (solid triangles) values 
of permittivity were considered. Error bars indicate the width 
of gap. vq does not depend on e for |e| > 10 4 . 
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FIG. 4. (a) Transmission through arrays of SRRs (dashed 
line) and through a LHM (solid line) with a complex per- 
mittivity for the metal, (b) Absorption versus frequency for 
the same system as shown in (a). Permittivity of metal is 
e = —3000 + i 100. Insert shows absorption in a compos- 
ite system with a frequency dependent e c ff and fi c s given by 
Eqns. ^ and (|). 
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FIG. 5. Phase of reflection for both SRR and LHM. 
e m = —1000 + i 18000. The vertical dotted line indicates 
the position of the center of the resonance gap vq. 
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FIG. 6. Transmission of EM waves through SRRs 
(dot-dashed lines), and LHMs (dashed lines). In (a) one SSR 
has been removed and its transmission is shown by a solid 
line with no appreciable difference from the periodic case. In 
(b) one metallic wire has been removed and its transmission 
is shown as a solid line. In this case there is a sharp dip in the 
transmission of the LHM with an impurity. In all the cases 
e m = —3000. For comparison, we also show as an inset in Fig. 
6b the transmission of waves through the lattice of wires only. 
A sharp resonance peak with a frequency Resonance ~ 13.65 
GHz appears when one wire is removed. Nine unit cells in the 
propagating direction are considered and the vacancies were 
created in the middle of the structure. 



